Arbitrary Control of Entanglement between Two Superconducting Resonators 
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We present a method to synthesize an arbitrary quantum state of two superconducting resonators. 
This state-synthesis algorithm utilizes a coherent interaction of each resonator with a tunable ar- 
tificial atom to create entangled quantum superpositions of photon number (Fock) states in the 
resonators. We theoretically analyze this approach, showing that it can efficiently synthesize NOON 
states, with large photon numbers, using existing technology. 
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The development of quantum coherent systems for in- 
formation processing has traditionally focused on quan- 
tum bits (qubits), in which information is stored in two 
quantum states of a system. Over the past two decades 
many physical systems have been devised in which qubits 
can be addressed and manipulated, including atoms, 
ions, photons, and solid-state systems pQ. However, 
recent experiments have demonstrated that supercon- 
ducting resonators — harmonic oscillators with a theoret- 
ically infinite ladder of states — can also be addressed 
and manipulated for quantum state storage and transfer 
[2]. These resonators have excellent coherence proper- 
ties, and would provide a promising alternative approach 
to large-scale quantum information processing. Future 
progress requires a theoretical study of how to efficiently 
generate entanglement in coupled networks of resonators. 

Recent experiments have achieved arbitrary control 
of a single superconducting resonator. In particular, 
Hofheinz et al. used a superconducting phase qubit to 
synthesize an arbitrary state of a single resonator [3J. 
While previously Fock states (states of definite photon 
number n) with n up to 20 had been generated [4], 
here superposition states were created and analyzed us- 
ing Wigner tomography [3J for photon states with n < 6. 
These states were synthesized using an algorithm de- 
veloped by Law and Eberly [5], originally designed for 
atomic cavity-QED systems. An important theoretical 
question is whether there exists a corresponding algo- 
rithm for the synthesis of an arbitrary quantum state of 
two resonators (a and b), of the general form 

N a N b 

l*> = J2 X! <=».,»»>«> ®|«6>- (1) 
n a = n b =0 

Among these states are the maximally entangled N- 
photon states 

1 N 

l*) = 7^yEl fc ^- fc )> ( 2 ) 

k—0 



and the so-called NOON states 

|*) = i(|iV„,0} + |0,iV i ,)). (3) 

States with the latter form can be used to beat the stan- 
dard quantum limit of measurements of phase (or fre- 
quency) [fy , and both represent generalizations of the Bell 
states of two qubits or the highly nonclassical N-particle 
Greenberger-Horne-Zeilinger state. 

In this Letter we present a realistic solution to the 
state-synthesis problem for two superconducting res- 
onators. We consider a tunable superconducting qubit, 
such as the phase [7] or transmon [8] qubit, coupled to 
two resonators with different frequencies, as shown in 
Fig. [T] This is described by the Hamiltonian 

h = w ,(t)|i)<i| + i(n(t)|i)<o| + n*(t)|o)<i|) 

+uj a a i a + LUbtfb (4) 
+g a (<7+a + er-af) + g b (a+b + a-tf) , 

where is the creation operator for the resonator of 
frequency u> a) W is the creation operator for a resonator 
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FIG. 1: Schematic circuit for generating entanglement be- 
tween two superconducting resonators. Resonator A (blue) 
has a fundamental frequency oj a /2n, while resonator B (red) 
has frequency 6J&/27T. These are each capacitively coupled to 
a tunable qubit (gray) with frequency u q /2-K, with coupling 
strengths g a and gb- The qubit is controlled by an external 
circuit (black). The theoretical results described in the text 
require uj a < uj q < u)b- 
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of frequency u>b, f2(f) is a possibly complex microwave 
field (in the rotating wave approximation), and g a and 
gb are fixed coupling strengths between the qubit and 
the resonators. Control of this circuit is exercised by 
modifying the time dependent qubit frequency u> q (t) by 
a "shift" pulse and applying Rabi pulses resonant with 
the qubit. By performing a sequence of these pulses, 
quanta can be created and transferred between the qubit 
and the two resonators. 

Note that this resonator-qubit-resonator system com- 
plements the qubit-resonator-qubit systems first studied 
in the resonant regime (uj q = uj a ) at NIST [5] and the 
dispersive regime (\uj q — uj a \ » g a ) at Yale [5] (note also 
the recent experiments [TOllllj). There have been several 
theoretical studies of this type of system [T2] , all towards 
the goal of generating entanglement between mesoscopic 
resonators. Here we solve the general problem of syn- 
thesizing an arbitrary entangled state, an important step 
towards quantum information processing with harmonic 
oscillator modes. 

A state-synthesis algorithm must program a sequence 
of pulses to prepare and transfer Fock states into the de- 
sired superposition state. Previous studies of this prob- 
lem for entangled states of motion for a single trapped 
ion [13] have shown that the synthesis of a general state 
with n a < -/V max and n b < N max requires a number of ele- 
mentary steps of order A^ ax , proportional to the number 
of coefficients in the expansion of the state vector, and 
schemes that achieve this scaling have been identified. 
However, none of these schemes can be directly applied 
to the problem presented above. These schemes all use 
sideband transitions and most use special two-mode in- 
teractions [H| specific to ion traps. While there are side- 
band interactions for resonators dispersively coupled to 
a qubit |15j . these interactions will be much slower; the 
transfer of a single photon has an effective Rabi coupling 

ofn cS ~ g \n\ 2 /uj 2 q m- 

For definiteness, consider the Fock-state diagram of 
Fig. [2]ja). Each node represents the two states of the 
qubit with photon numbers (n a ,n(,), or the quantum 
state |g, n a , rif,} (where the qubit state is q = or 1). The 
resonant interaction of the qubit with each resonator, 
H a = g a { a + a + <J-a)) and Hi, = gb(<J+b + <r_6^), leads 
to horizontal and vertical transitions, respectively, illus- 
trated by the solid and dashed lines. These resonant 
interactions are fast, efficient, and provide nearest neigh- 
bor transitions in the Fock-state diagram. The addition 
and control of individual photons requires an indepen- 
dent state-selective qubit rotation. By addressing the 
qubit in between the resonant cavity interactions, it is 
possible to access the entire state space. 

To achieve selective manipulations of the quantum sys- 
tem, we use the photon-number-dependent Stark shift 
[T7j . For our system, this implies that a qubit operated 
in the dispersive regime will undergo Rabi oscillations 
from |0,n a ,nf,) —¥ |1, n a ,nt,) when the drive frequency 



satisfies 



-(2n a + l) + 



9l 



- Ub 



-(2n b + l). (5) 



We choose to set Aw = 2g 2 /(uj q — u> a ) = —2g 2 /(oj q — uif,). 
This can always be achieved for a qubit with a tunable 
frequency such as the phase or transmon qubit. This 
allows us to simplify our Rabi pulses to frequencies 



ui n = u) q + nAw, 



(6) 



which selects those states with n a — rib = n, where n is 
an integer (n — is shown in Fig. [5] (a)). Note that this 
choice requires uj a < io q < u>b and < <7^/(w g — u>b) 
(to avoid nonresonant transitions) . By choosing different 
values of u n one can address each of the "diagonals" of 
the Fock-state diagram. 

Direct simulations of Rabi oscillations using the full 
Hamiltonian verify this approach, as shown in Fig. [2jb) 
for n = 2. Each block represents the maximum of 
the transition probabilities \(l,n a ,ni,\'i'(t))\ 2 , calculated 
with the initial condition \^(t = 0)) = |0,n o ,n(,). The 
transition probabilities are large along the diagonal, de- 
creasing significantly for neighboring Fock states. A sim- 
ilar effect known as "number splitting" was experimen- 
tally demonstrated for a qubit coupled to a single res- 
onator |18) . and recently used in an experiment with a 
qubit coupled to two resonators [llj . 

We now show how these three interactions can be used 
to achieve the synthesis of arbitrary joint states of two 
resonators. This is accomplished by the following se- 
quence of operations 



U 




(7) 



with 



N a N a 

I 1 J AjR^ and U bJ = \[ B jk R bdk . (8) 



k=0 



Here we have defined Aj 



exp(-iH a t aJ ), Bj 



cxp(— i-fffci{,jfe), and the single-qubit rotations R a j and 
Rb,jk- This sequence can be physically realized by shift- 
ing the qubit into and out of resonance with resonators 
a and b, interleaved by the Stark-shifted qubit rotations 
described above. The parameters (t a> j, tb,jk, Ra,j, and 
Rbjk) m this operation must be chosen to satisfy 

N a N b 

|*) = C/|0,0,0) = |0) ® 2 H c„ a ,„> a ,n b ), (9) 

n a = l n& = l 

where c na>ni are arbitrary coefficients. To determine the 
precise sequence of operations for a given two-resonator 
state 1^), one solves for the inverse evolution: 
l 

u^\*) = ul J] [/^.l*) = |0,0,0). (fO) 

j=N b 
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FIG. 2: Schematic set of operations to generate an arbitrary state of two resonators. In this Fock-state diagram, the state 
\q, n a ,nb) is represented by the node at location (n a ,n\,). (a) Interactions lead to couplings between these states, indicated by the 
arrows. Three key interactions are used: A transfers quanta between the qubit and resonator a (solid lines), B transfers quanta 
between the qubit and resonator b (dashed lines), and R (curved arrows) rotates the qubit for states with n a — rib = n (here 
with n — 0, see text), (b) Numerical simulations of Stark-shifted Rabi oscillations for uj a /{2-n) — 6.3 GHz, W(,/(2tt) = 7.7 GHz, 
co q /(2ir) = 7 GHz, g a /{2Tr) = gb/(2-ir) = 70 MHz. The Rabi oscillations are driven at cj d /(27r) = 7.025 GHz and Q/(2n) = 7 MHz 
(with n a — ni, — 2). Each block corresponds to the maximum probability of the transition \Q,n a ,rib) — > |1, n a ,rib) (see text). 



Each step of this inverse sequence can be mapped onto 
a state transfer in the Fock-state diagram. In order to 
solve this problem, one must show that photons can be 
consistently removed from the system. Our approach ac- 
complishes this in the following way. The sequence of 
U% a unitaries moves the system along the vertical paths 
(B) in the Fock-state diagram. Each step removes a pho- 
ton from row rib = j of the Fock-state diagram. This is 
done by marching (from right to left) along the columns 
n a = k, with Bjk transferring the amplitude in row j 
to row j — 1, after which the Stark-shifted single-qubit 
operations Rb,jk (with frequencies ui q + (fc — j + l)Au) 
rotates the amplitude to state |0, k,j — l). After all of the 
photons have been removed from the columns in row j, 
the sequence repeats for row j — 1. Each time through, 
population is transferred towards rib = 0. Once there, 
the Ul sequence moves population along the horizontal 
paths (A) to n a — nt — 0, or |0, 0, 0), thus solving Eq. 
(10 1. This completes the algorithm. 



total time is approximately given by 



The total number of steps matches the optimal effi- 
ciency of the ion-trap proposals discussed above, but here 
using resonant interactions and the Stark-shifted Rabi 
pulses. Each step involves the rotation of an effective 
two-state system whose amplitudes are known (from the 
original c nat n h ), as in the original Law-Eberly scheme [5]. 
By counting the number of operations in U, we find that 
the general sequence requires N a A unitaries, (N a + l)Nb 
B unitaries, and N a + (N a + l)Nb Rabi pulses. Thus, the 



(JV +1)(JV 6 +1) 



SI 



N a 



<N a +l) 



N b 

j- 



(11) 

Note, however, that for states such as the NOON state 
we can achieve an even greater efficiency. For these states 
one need not transfer amplitude over the whole diagram, 
but only along certain paths, leading to a sequence with 
only a linear number of steps. Consider the sequence of 
operations shown in Fig. [3| whose steps are detailed in 
Table |lj This sequence requires a linear number of oper- 
ations as opposed to the quadratic scaling of the general 



procedure described above. In fact, we find that Eq. ( 11 ) 
can be reduced to 



NOON 



N n + N h 



7T 
ft 



(12) 

We now estimate the time required to generate a 
NOON state. We consider a qubit with uj q j (2tt) = 
6.5 GHz and resonators with u a /(2ir) — 6 GHz, 
Wb/(27r) = 7 GHz and a coupling of g a /{2Tr) = gb/(2ir) = 
150 MHz (similar to recent experiments [19 J. For the 
state-selective Rabi pulses, numerical simulations (not 
shown) show that Q/(2?t) = j^g 2 a /{io q - u a ) = 22 MHz 
produces an error of a few percent. Using Eq. (12 1 we 
estimate that NOON state generation with N a = Nb = 8 
will take only 360 ns. Smaller couplings (as in Fig. 2) can 
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FIG. 3: Algorithm to generate the state |*> = |3, 0) + |0, 3) of 
two coupled resonators. The sequence of operations is detailed 
in Table 1. The horizontal, vertical, and curved transitions 
are the interactions A, B, and 7? (see text). 



still achieve N a = Nf, = 3 in 410 ns. These times com- 
pare quite favorably to the coherence times of both qubits 
and resonators, which are now consistently greater than 
0.5 fis [3111 [TH]. Most of the time is for the high-fidelity 
state-selective Rabi pulses. Faster rotations should be 
possible by using specially shaped pulses [5D]. Recent 
experiments [21] using such pulses show that quantum 
algorithms are ultimately limited by the coherence times. 

Other experimental issues may arise in this procedure. 
First, there will be modifications to the rotating wave 
and dispersive approximations (used to derive Eq.([5])) 
for large couplings and photon numbers. These can be 
addressed through pulse shaping or optimal control the- 
ory approaches. Second, we have ignored the dynamical 
phases that arise when the qubit is shifted between fre- 
quencies. These phases can be corrected by including 
brief pauses between the Rabi and shift pulses [3J IT9] . 
and will not significantly affect the overall time needed for 
state preparation. A full simulation including these ef- 
fects, including decohcrence, will be performed elsewhere, 
but the estimates given above are quite promising. Fi- 
nally, verifying the two-resonator state may require addi- 
tional qubits for readout using Wigner tomography [3J to 
probe the coherence of the two resonators. However, effi- 
ciently manipulating and measuring entangled resonators 
appears experimentally possible. 

In summary, we have presented a method to synthe- 
size an arbitrary quantum state of two superconduct- 
ing resonators. This method combines state-selective 
Rabi oscillations (using Stark shifts) with linear cou- 
plings of each resonator to a tunable artificial atom. We 
have shown that this approach can efficiently synthesize 
NOON states, with large photon numbers, using existing 
technology. More generally, this approach can be applied 
to many types of coupled qubit-resonator systems, and 



TABLE I: NOON State-synthesis procedure 



Step Parameters Quantum State 



Ra,l 


Qtqa, 1 


= Tv/2, LU d = LO 


1 0,0,0) - i 


1,0,0) 


Ai 


9ata,l 


= tt/2 


|0,0,0)- 


0,1,0) 


Ra,2 


Qtqa,2 




1 0,0,0) +i 


1,1,0) 


M 


g a ta,2 


= tt/(2V§) 


jo, 0,0) + 


0,2,0) 


Ra,3 


Qtqa,3 


= 7T, Ud = L02 


|0,0,0)-i 


1,2,0) 


A 3 


g a ta,3 


= tt/(2v / 3) 


10,0, 0> - 


0,3,0) 


R b ,i 


Qtqb,l 


= 7T, UJd = Wo 


-i|l,0,0) - 


0,3,0) 


B x 


gttb,i 


= tv/2 


-|o,o,i)- 


0,3,0) 


Rb,2 


Qtqb,2 


— TV, UJd = W-l 


i|l,0,l)- 


0,3,0) 


B 2 


gttb,2 


= V(2\/2) 


|0,0,2)- 


0,3,0) 


Rb,3 


Qtqb,3 


= 7T, UJd = UJ-2 


— 0, 2) - 


0,3,0) 


B z 


gbtb,3 


= 7r/(2\/3) 


-|0,0,3) - 


0,3,0) 



opens up an important path towards quantum informa- 
tion processing with superconducting oscillators. 
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